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ABSTRACT 

Aims. We study the expected properties of starburst galaxies in order to provide the point of reference for interpretation of high-z 
galaxy surveys and of very metal-poor galaxies. We concentrate mainly on the UV characteristics such as the ionizing spectra, the 
UV continuum, the Lyo- and He n /tl640 line and two-photon continuum emission. 

Methods. We use evolutionary synthesis models covering metallicities from Pop III to solar and a wide range of IMFs. We also 
combine the synthetic SEDs with the CLOUDY photoionization code for more accurate predictions of nebular emission, and to study 
possible departures from case B assumed in the synthesis models. 

Results. The ionizing fluxes, UV continuum properties, and predicted Lya- and He n /11640 line strengths are presented for synthesis 
models covering a wider range of parameter space than our earlier calculations. Strong departures from case B predictions are obtained 
for Lya and 2y continuum at low metallicities. At low nebular densities both are shown to be enhanced proportionally to the mean 
energy carried by the Lyman continuum photons emitted by the ionizing source. Larger Lyo- equivalent widths are therefore predicted 
at low metallicity. The He n /ll 640 line can be weaker than case B predicts (in terms of flux as well as the equivalent width) due to its 
ionization parameter dependence and to the enhanced underlying 2y continuum. 

Conclusions. Our results have implications for the interpretation of star-forming metal-poor and/or high redshift galaxies, for galaxies 
among the Lyo- emitters (LAE) and Lyman Break galaxy (LBG) populations, and for searches of Population III stars in the distant 
Universe. 
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lutionary synthesis models and has demonstrated the impor- standard evolutionary synthesis are. 
tance of nebular emission (lines and contin ua) at low metallic- 
ity. These and o ther st udies (see e.g iTumlinson & Shulll l2000: In deed, a shortcoming of evolutionary synthesis models such 



Tumlinso nttall l200ll Bromm et al. 2001^ have in particular as th e ones mentioned above is that they calculate nebular emis 

sion in an approximate manner assuming simplifi ed physics. 
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2006), and constant emissivities for adopted constant values of 
the electron temperature and density in the H n region sur- 
rounding the starburst. In fact, as demonstrated in this paper, 
significant departures from case B are expected at low metal- 
licities leading to stronger Lya emission, and the strength of 
nebular He n emission predicted by full photoionization mod- 
els can be reduced with respect to simple recombination the- 
ory. Indeed these physical effects, related to an increased im- 
portance of collisional effects at low metallicity due to lower 
radiative cooling and harder ionizing spectra - for Lya - and 
due to competition between H and He for ionizing photons - for 
the intensity of He n/H - have been known for a while in stud- 
ies of metal-poor H n regions (cf. iDavidson & KinmanH l985: 
Stasiriska & Schaerer 1999; Luridiana et al an d plane tary 

nebulae dStasiriska & Tvlendalll986l) . lPanagial d2002T 120051) has 
recently explored photoionization models for primordial nebu- 
lae. However, the importance of the above effects for the UV 
emission lines has so far not been thoroughly examined, in par- 
ticular in the context of metal-poor and distant starburst galaxies 
and using up-to-date evolutionary synthesis models. The pho- 
toionization models presented here, combined with our evolu- 
tionary synthesis models, are intended to provide a framework 
within which to improve our knowledge of primeval galaxies 

and related objects. 

Another limitation of the synthesis models of ISchaererl 
(120021120031) concerns the initial mass function (IMF). For sim- 
plicity, three different choices of the IMF were adopted for the 
bulk of the calculations in these papers. However, different IMFs 
have been suggested in other studies related to Pop III and early 
stellar generations, and considerable uncertainties remain on the 
true shape of the IMF in the early Universe and its dependence 
(or not) on physical parameters. To enable the examination of 
the effects of a broader choice of IMFs on the expected observ- 
able properties o f starbursts, we here extend the calculations of 
Schaerei] (120031) to eight different IMFs. The resulting model 
grids, available in electronic format, should provide state-of-the- 
art predictions for the interpretation of high redshift galaxies, to 
estimate their contribution to cosmic reionization, and for other 
topics. 

The paper is structured as follows. In Sect. [2] we describe 
the input physics and the model calculations with our evolution- 
ary synthesis code and with the photoionization code CLOUDY. 
The predictions from the synthesis models concerning the UV 
continuum, the ionizing flux, Lya and He n A 1640 emission are 
presented in Sect. [3] In Sect. [4] we discuss the results from the 
photoionization models using black body spectra, explain the de- 
viations from case B and provide simple formulae to describe 
these effects on Lya. In Sect. [5] we show how to connect realis- 
tic SEDs with results from photoionization models using black 
body ionizing spectra. Our results and several implications are 
discussed in Sect. [6] The main results are summarised in Sect. [7] 



2. Modeling techniques 

2.1. Synthesis models 

We have used the evolutionary synthesis code of 
ISchaerer & Vaccal dl998l) . with the physical ingredients 
(stellar tracks, atmospheres, an d prescrip t ions f or nebular line 
and continuum emission) from Schaerer (2003). In particular 
these models allow us to predict the integrated properties of 
stellar populations at all metallicities from zero (Population III) 
to "normal", solar-like metallicity. The computation s have been 
done for the same metallicities as in ISchaererl d2003h . 



2.1 .1 . Stellar initial mass function 

The ma in extension pre sented here with respect to the calcula- 
tions of ISchaererl d2003) concerns different assumptions regard- 
ing the stellar IMF. A wide range of IMFs has bee n cons i dered , 
including power-law IMFs, such as the Salpeter or rScalol (fl986) 
IMFs, log-normal IMFs, and the iLarsonl d!998l) IMF. The cor- 
responding parameters are summarised in Table Q] The stellar 
mass range is defined by the lower and upper mass cut-offs, Mi ow 
and M up respectively, a being the slope of the power-law. The 
log-normal IMFs are described by the characteristic mass M c 
and its dispersion <x. The cases computed here corres pond to the 
cases A, B, and E in the chemical evolution study of Tumlinson 
(2006). The lLarsonl dl998) IMF is described by a single parame- 
ter, its characteristic mass M c . We h ave computed one such case, 
assuming the same value of M c as ICiardi et all d200l in their 
reionization calculations. Note, that at Z > 0.0004 = 1 /5OZ the 
upper mass cut-off is set to M up = 100 or 12OM for all IMFs, 
since tracks for more massive stars are not available. The quan- 
tities discussed here are insensitive to assumptions on the IMF 
at low masses. Our absolute quantities may therefore simply be 
rescaled to other IMFs including e.g. an extension below 1 M . 

Current knowledge suggests that the IMF for massive stars 
is close to Salpeter with an upper limit of M up ~ 100-120 M 
from solar metallicity down at least to ~ 1/50 Z , and that a qual- 
itative shift of the IMF towards higher characteristic masses oc- 
curs below a critical metallic ity of the order of Z cr i t ~ 1O~ 5±1 Z 
(ISchneider et al.ll2002ll2003l) . 



2.1.2. Star formation histories 

For each metallicity Z and IMF we have computed evolutionary 
synthesis models for the two limiting cases of 1) an instanta- 
neous burst, and 2) constant star formation (CSFR). Results for 
other star formation histories can be derived from the electronic 
files for the simple stellar population (burst) models (see Sect. 
12.1.4b . In both cases the calculations have been carried out with 
a small time step (0. 1 Myr) to ensure the accuracy of the time in- 
tegrated quantities for the CSFR case. The calculations have to 
be carried out up to ages of 1 Gyr. This covers the allowed ages 
and of galaxies at redshifts z ^ 5.8, of interest here, as well as 
sufficiently long star formation timescales to reach equilibrium 
in various observable properties (cf. below). 

2.1.3. Nebular emission 

To include nebular emission (recombination lines and continuum 
processes) in our synthesis models we initiall y make the fol- 
lowin g "standard" simplifying assumptions (see ISchaererf 2002, 
2003): ionization bounded nebula, constant electron tempera- 
ture and density (T e , n e ), and case B. Case B in particular as- 
sumes that the recombinations to the ground-state immediately 
yield locally another ionization, and that photoionizations oc- 
cur only from the ground-state. As we will see below, the latter 
may not be true in very metal-poor nebulae, leading to signifi- 
cant changes in the predicted spectrum of hydrogen. With these 
assumptions both the H and He recombination lines as well as 
nebular continuum emission (including free-free and free-bound 
emission from H, neutral He, and singly ionized He, and two- 
photon emission of H) are fully specified and their luminosity is 
proportional to the ionizing photon flux Q in the appropriate en- 
ergy range. To reflect to first order the changes of the conditions 
in the H n regions with metallicity, the value of the line lumi- 
nosity coefficient and nebular continuous emission coefficients 
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Table 1. Summary of IMF model parameters . Mode l ID is the label used in the Figures, colour code the colour and linestyle. Note: 
the definition of M c and <x is as in Tumlinson (2006). In particular sigma is the variance in ln(m), not log(m)! 

Model ID colour code Mi ow M up a M c a reference/comment ID in files 



Salpeter 


black 


1 


100 


2.35 






A in Schaerer (2003) 


S 


B 


green, dashed 


1 


500° 


2.35 






B in Schaerer (2003) 


B 


C 


cyan, dashed 


50 


500° 


2.35 






C in Schaerer (2003) 


C 


Scalo 


blue 


1 


500° 


2.7* 






Scalo (1986) 


Sc 


TA 


red 


1 


500° 




10. 


1.0 


A in Tumlinson (2006) 


TA 


TB 


magenta 


1 


500° 




15. 


0.3 


B in Tumlinson (2006) 


TB 


TE 


yellow 


1 


500° 




60. 


1.0 


E in Tumlinson (2006) 


TR 


L05 


blue, dashed 


1 


100 




5. 




Larson (1998) 


10 



" For metallicities Z > 0.0004 = l/50Z o , we adopt M up = 120M o , the maximum mass for which Geneva stellar evolution tracks are available. 
* Power-law exponent for M > 2M Q . 



at metallicities Z/Z < 10 are changed as in ISchaererl (120031) . 
More precisely we adopt T e = 30 (20) kK for lines (continua) 
Z/Z < 10~ 3 and T f =10 kK for higher metallicities, and a low 
ISM density n e = 100 crrT 3 . 

For the Lya luminosity we have, with the assumptions just 
spelled out: 



L B (Lya) = (1 - / esc )G(H) x hv Lya x 



T eff 



(1) 



where the index 'B' stands for case B, a B = ar^L + i s me 
total case B recombination coefficient, and / esc is the escape frac- 
tion of ionizing photons out of the H n region (or galaxy). In all 
Figs, shown in this paper we assume / esc = 0. In typical con- 
ditions a^p/ag ~ 0.6-0.7. In other words approximately 2/3 of 
the Lyman continuum photons give rise to the emission of a Lya 
photon, the assumption commonly made. Similar relations also 
hold for other recombination lines, such as He n A 1640, which 
is of special interest here. In our "standard" synthesis models we 
simply assume 



L B (Lya) = Q(H) Xc u 
L B (He IU1640) = g(He + ) x c 2 , 



(2) 
(3) 



with ci = 1.04 x 10~ n erg, and c 2 = 5.67 x 10~ 12 (6.04 x 10~ 12 ) 
erg for Z < 1/50 Z (> 1/50 Z ), and for / esc = 0. The atomic 
data is from H ummer & Storey (1995) for low densities (cf. 
Schaerer 2003), and g(H) and g(He + ) are the ionizing photon 
flux (in photon s ) above 13.6 and 54 eV respectively. 

Continuous nebular emission including free-free and bound- 
free emission by H, neutral He, He + , and He +2 , as well as the 
two-photon con tinuum of hydrogen is included as described in 
ISchaererl (120021) . assuming T e = 20 kK for Z/Z < 5 . x 1 4 and 
T e = 10 kK otherwise. 

As we will show below (Sect. @), a proper treatment of all 
relevant processes leads to significant deviations from case B at 
very low metallicities, increasing in particular the Lya luminos- 
ity, L(Lya). In this case the L(Lya) can be rewritten as 



L(Lya) = L B (Lya) x P x — , 



(4) 



where P and / co u are terms describing the mean photon energy in 
the Lyman continuum, and accounting for for collisional effects 
at high density. To compu t e the Lya e quivalent width W(Lya) 
we proceed as in ISchaererl d2002l 12003). where we use the con- 
tinuum flux at 1215.67 A obtained from linear interpolation of 
the total (stellar + nebular) continuum (in log) between 1190 



and 1240 A, chosen to avoid underlying stellar Lya absorp- 
tion and other absorption lines if present. The predicted stellar 
Lya absorption is small compared to t he emission, except for 
"post -starburst" phases (see e.g. Fig3. in ISchaerer"&~ Verhamme 
2008). 



2.1.4. Model output 

Our evolutionary synthesis code predicts a large variety of ob- 
servable and related quantities derived from the detailed syn- 
thetic spectra. Here we focus on mostly on quantities describ- 
ing the spectrum in the Lyman continuum, the UV (rest-frame) 
spectrum, as well as the Lya and He n A164Q emission lines. 
The full set of model results, including also numerous other 
quantities not discussed in this paper, are available in elec- 
tronic format upon request to one of the authors (DS), on 
http : //obswww . unige . ch/sf r or via the CDS. 

2.2. Photoionization models 

To predict more accurately the nebular emission from starbursts 
and to investigate possible departures from the simplified as- 
sumptions made in our synthesis m odels we use the p hotoioniza- 
tion code CLOUDY version 08 dFerland et al.lll998l) . The mod- 
els we consider are ionization-bounded with a closed, spherical 
geometry and constant density. These assumptions imply in par- 
ticular that all ionizing photons are absorbed in the H n region, 
i.e. /esc = 0. In certain circumstances, especially in high red- 
shift galaxies, a fraction o f the Lyman continuum photons are ex- 
pected to escape (see e.g.lGnedin et al.ll2008t IWise & Cenl l2009: 
iRazoumov & S ommer-Larsenl l2009l )7To first order the results 
obtained in this paper can simply be rescaled to such cases, as 
discussed below. The main input parameters of our models are: 
the spectral energy distribution (SED) of the ionizing source, the 
nebular density (given by n^), the hydrogen number density), 
the ionization parameter (£/), and the nebular metallicity (Z ne b). 
For the SED we adopt black-body spectra described by T/,t and 
SEDs from our evolutionary synthesis models. The ionization 
parameter (at the inner edge of the cloud) is defined as: 



U 



Q(H) 



Anr 1 x n(H) x c 



(5) 



where r m is the inner radius of the nebula which has been kept 
constant in our models (10 17 cm) and c is the speed of light. The 
small inner radius used in our calculations results in a sphere-like 
(not shell-like) geometry of the nebulae. The ionization param- 
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SFR = const @ 1 Gyr ZAMS values 



13.5 



12.5 




-4 -2 
log(Z/Z a ) 



Fig. 1. Relative output of hydrogen ionizing photons to UV continuum light, measured at 1500 A restframe, Qh /£i500> as a function 
of metallicity for constant star formation over 1 Gyr (left panel) and very young bursts (right panel). 2///Li5oo is given in L,\ units 
on the left side of each panel, and in L v units on the right. Results for different IMFs are shown using the colour codes and labels 
summarised in TableQ] The shaded area indicates the critical metallicity range where the IMF is expected to change from a "normal" 
Salpeter-like regime to a more massive IMF (see text). 



eter can change throughout the nebula (decrease towards outer 
parts). 

First, we have calculated the grid of photoionization models 
using the photoionization code CLOUDY covering log(£/) = -4, 
-3, -2, -1; T hb = 40,000-150,000 K; log(n(H)) = 1, 2, 3, 4 cirr 3 
for a primordial nebula and a number of higher metallicities. In 
total 192 models were computed for each metallicity. In order to 
keep the same ionization parameter while changing the density 
of the gas, we adapt the number of ionizing photons (keeping 
the shape of the SED). The metallicity is defined by scaling the 
solar abundance pattern. 

Note that the electron temperature is neither constant nor 
fixed in our models. Its spatial distribution results from the com- 
putation of each photoionization model and it is a function of 
depth in the nebula. In practice it depends on all the parameters 
that are being investigated (Tt,i„ tin, U, Z„ e b). For the coolest 
primordial model (the coolest black body, the lowest density and 
the lowest ionization parameter) the temperature in the inner part 
of the cloud is around 12,000 K and for the hottest ones it reaches 
~ 38,000 K. 

Selected models were subsequently computed using the 
SEDs from the synthesis models described above. 

3. Predicted UV properties from synthesis models 

We now present and discuss one-by-one the main predictions of 
our synthesis models for different IMFs (see Table [TJ, for metal- 
licities from zero (Pop III) to solar, and for two different limiting 
star-formation histories (bursts and SFR=const). Since proper- 
ties of stars below Z S> 10~ 9 (i.e. Z/Z ^ 10~ 73 ) essentially 
converge to those of meta l-free star s we a ssign this metallicity 
value to Pop III stars, as in lSchaereri (120031) . 

Note that all UV continuum predictions from the synthesis 
models described in this Section are based on the simplified as- 



sumptions spelled out above to compute nebular emission. This 
implies in particular that at low metallicities the contribution 
from the two-photon continuum process should be higher, in- 
creasing thus e.g. the predicted UV luminosity at 1500 A, for 
the reasons discussed in Sect. [4] 

3. 1 . Ionizing photon production 

A quantity of interest, e.g. to determine the contribution of galax- 
ies to cosmic reionization, is the the relative output of hydro- 
gen ionizing photons to observable UV light. Here we provide 
Qh/L\5oo, where the Lyman continuum flux Qh is expressed in 
units of photon s _1 , and the UV restframe luminosity at 1500 
A is Z4500 in Lx (erg s _1 A -1 ) or L v (erg s _1 Hz -1 ) units Q 
Alternatively, to derive the Lyman continuum output per unit 
stellar mass (or per baryon), Qh/^isoo discussed here can be 
combined with the ratio SFR/Z4500 given below, or can directly 
be derived from the available data files. 

In Fig. Q] we show Qh/^-isoo as a function of metallicity 
for constant star formation over 1 Gyr (CSFR, left panel), and 
very young {% Q-A Myr) populations (right panel). As expected 
Qh/L\5oo increases with decreasing metallicity, since the ioniz- 
ing flux depends very strongly on the effective stellar tempera- 
ture and hence increases more rapidly than the UV luminosity. 
The IMF dependence also behaves as expected, with the IMFs 
favouring the most massive stars showing also the highest the 
Qh/Lisoo ratios, since Qh increases more rapidly with stellar 
mass than the UV luminosity. Notable is actually the increase 
of Qh/Lisqo from Z to ~ 1O~ 4 Z , where no major change of 
the IMF is expected (and a Salpeter IMF is favoured). For CSFR 
and for a fixed IMF, the increase of the relative ionizing power 

1 The transformation is log(Q H /L A (1500A)) = log(e„/L v (1500A))- 
12.12. 
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time evolution for bursts/CSFR models and different IMFs 
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Fig. 2. Temporal evolution of (2///^i500 for selected IMFs 
(Salpeter, TB, TE, colour-coded as in previous Figs.) and metal- 
licities (Pop III, Z , labeled). The solid curves show the time 
evolution for constant SFR models towards their equilibrium 
value, the dashed curves instantaneous burst models. See text 
for discussion. 

SFR = const @ 1 Gyr (thick), ZAMS (thin) 
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Fig. 3. Mean ionizing photon energy in units of 13.6 eV as a 
function of metallicity, shown for SFR=const (thick lines) and 
for the ZAMS (thin lines). Results for different IMFs are shown 
using the same colour codes as in Fig.Q~](cf. Table[T|). Lined and 
labels on the left indicate the corresponding blackbody temper- 
atures. Z co u shows the approximate metallicity limit below with 
collisional effects lead to significant departures from case B. 



from solar metallicity to Pop III is typically ~ 0.4-0.5 dex, i.e. a 



factor 2 to 3. When considering an IMF change from Salpeter to 
a massive IMF (i.e. all cases except Salpeter and Scalo) the in- 
crease of Qh/Lxsm is larger, approximately 0.6 to 1 dex between 
solar and zero metallicity. Only for the "TB" IMF, a narrow, log- 
normal mass function peaked at M — 15M Q , we find a more 
extreme dependence on metallicity. This is precisely due to the 
fact that this IMF singles out a narrow mass range, instead of 
averaging the metallicity dependence of stellar properties over a 
larger interval in mass. 

The right panel of Fig.[T]shows a narrower range of Qff/^isoo 
for zero age or very young {% Q-A Myr) populations. This is nat- 
ural, since in this case no "average" is made over populations 
of very different stellar ages and hence over strong variations of 
stellar parameters. More important is the fact that higher values 
of Qh I 'Z4500 are obtained for young populations. This is mostly 
due to the fact that such populations emit a lower UV luminos- 
ity per unit SFR since a longer timescale is needed to reach the 
"equilibrium value" of the UV output (cf. below). The ZAMS 
values shown here correspond to the maximum of Qh/^isoo ex- 
pected for stellar populations of different ages and SF histories. 

To illustrate this dependence on the SF timescale, 
Qh / Lisoo(t) is shown in Fig.|2]for selected IMFs and metallici- 
ties. These curves show the smooth transition from the predicted 
"ZAMS" to the CSFR values over timescales from ~ 10 7 yr for 
massive IMFs (e.g. TE, TB) to ~ 0.4-1 Gyr for the Salpeter IMF. 
Note also that the timescale for UV properties to reach equilib- 
rium increases with decreasing metallicity, due to the higher ef- 
fective temperatures on the ZAMS at low Z. In short, we caution 



SFR=const @ lGyr 




log(Z/Z e ) 

Fig. 4. Dependence of the SFR(UV) calibration on metallicity 
and IMF. Shown is the SFR per unit UV luminosity at 1500 
A in units of M Q yr~' per (erg s A -1 ) on the left y-axis, or 
in M yr per (erg s Hz ') on the right y-axis. Same sym- 
bols and colours as in previous Figures (cf. Table [TJ. These val- 
ues assume constant SF over 1 Gyr. The dotted horizontal line 
show the value of the Kennicutt (1998) SFR(UV) calibration 
rescaled to a Salpeter IMF with Mi ow = 1M for comparison. 
Above/below this line the SFR deduced from the UV luminosity 
is higher/lower. 
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that the relative ionizing photon to UV ratio Qh/L\5Q0 may be 
uncertain by a factor of ~ 4 depending on the SF timescale (for 
constant SF), or even more for bursts. 

Finally, it should be noted that at low metallicity the con- 
tribution of the two-photon continuum may be larger, as shown 
later in Sect. [4] leading to somewhat lower values of Q/f/^isoo- 
For constant SF this amounts to a decrease of (2///£i500 by 
~ 0.2 dex for Pop III and the most extreme IMFs (TE, C), and 
smaller changes otherwise. For zero metallicity populations on 
the ZAMS Q H /Li 500 should be reduced by ~ 0.15-0.3 dex for 
all IMFs, and less at higher metallicity. 

3.2. Properties of the ionizing spectra 

The properties of the ionizing spectra, such as their har dness, de- 
tailed shape and others have already been discussed in lSchaererl 
(2003) and shall not be repeated here. For example, the hard- 
ness g(He + )/g(H), measured by the ratio of He + (> 54 eV) to 
hydrogen ionizing (> 13.6 eV) photons, we predict from our 
new models are already bracketed by the the values predicted in 
ISchaererl d2003l) (see their Fig. 5) for the Salpeter and the "C" 
IMF. 

An interesting quantity describing the ionizing spectra is the 
average energy of the photons emitted in the Lyman continuum 
E(Q(HJ) (see definition in Eq. |6j. This quantity and its depen- 
dence on metallicity and IMF is plotted in Fig.[3]for constant SF 
(thick lines) and for the ZAMS (thin lines). The corresponding 
blackbody temperatures 7bb with the same mean ionizing pho- 
ton energy are also shown for illustration. Typically E(Q(H)) 
is found to ~ 1.5-2.5 times 13.6 eV, the ionizing potential of 
neutral hydrogen, and its behaviour with IMF, metallicity, and 
age behaves as expected. These values correspond to a range of 
blackbody temperatures from ~ 50 to 120 kK for the hardest 
spectra. 

On this Figure we also indicate the approximate metallicity 
limit Z co ii, below which collisional effects lead to significant de- 
partures from case B, as shown below. In this metallicity range 
E{Q(H)) can also be used to compute more accurately the intrin- 
sic Lya emission line strength (see Eq.[8]l. 

3.3. SFR calibrations from the UV continuum 

Figure |4]illustrates the variation of the UV luminosity for CSFR 
as a function of metallicity and for the different IMFs. Plotted 
is the conversion factor c, defined by SFR= c x L v , where L Y 
is the UV luminosity in units of erg s _1 Hz -1 , and SFR is the 
star formation rate in M yr _1 . As expected our model with 
Salpeter IMF agrees well with the widely used calibration from 
Kennicutt (1998) at Z after rescaling the latter by a factor 2.55 
to account for our adopted value for the lower mass cut-off 
(Mi ow = 1M 2- This Fig. clearly shows that in most cases the 
use of the Kennicutt calibration at low metallicity may overesti- 
mate the SFR, given the higher intrinsic UV output of such stel- 
lar populations. However, this may not be realistic since it relies 
on the assumption of CSFR over a long timescale (S 10 8A 9 yr). 
In younger populations the UV luminosity per unit SFR is lower, 
and hence c and the determined SFR higher (e.g. Schaerer 2000). 

In fact, the non-monotonous behaviour of SFR/Z4500 with 
metallicity observed in Fig. [4] is due to the dependence of the 



2 More explicitely the Fig. shows logc, where log SFR = logc + 
logLisoo + log Cm, and c M = 2.55 for the IMF adopted by Kennicutt 
(Mow = 0.1 Ms), or c M = 1 for Mi ow = 1 M Q . 



stellar contribution to the total UV output at this wavelength. 
Indeed, the stellar UV output increases with decreasing Z down 
to ~ 1/50 Z/Z , due to the decrease of the average temper- 
ature of stars over their lifetime. At even lower metallicities, 
however, their UV output (per unit SFR) decreases since the 
bulk of their flux is emitted at <k 1500 A(cf. Fig. 2 in lSchaererl 
2003). This implies, for a fixed IMF, a re-increase of SFR/Z4500 
at Z/Z ^ 1 /50, which is only somewhat moderated by the in- 
creasing nebular contribution. Indeed, the latter contributes typ- 
ically ~ 10-40% of the flux at 1500 A at SF equilibrium (see 
Fig. [8) In other words neglecting the nebular continuum would 
lead to differences of ~ 0.05-0.15 dex in the SFR calibrations. 
At low metallicity the contribution of the two-photon continuum 
may be larger, as shown in Sect. |4] leading to somewhat lower 
values of SFR/L1500 than shown in Fig. [4] For Pop III and the 
most extreme IMFs (TE, C) this implies a downward revision of 
~ 0.2 dex. 

3.4. Predicted Lya emission 

The Lya equivalent widths predicted by our standard models 
(using Eq. [2]i for all IMFs and metallicities are shown with thin 
lines in Fig.|5]for constant star-formation (left panel) and for the 
ZAMS (right), the latter representing the maximum of W(Lya) 
for each IMF. The dependences of W(Lya) are as expected from 
Qff/Lisno and the values shown here bracket those already pre- 
sented in ISchaererl d2003l) (with more limited variations of the 
IMF). To illustrate the departure from case B found at low metal- 
licity from photoionization models (see Sect. |H, we also show 
W(Lya) computed from Eq.|4]in the low density regime (i.e. for 
^oii = 2/3, thick lines). Here, the equivalent widths are a fac- 
tor ~ 1.5-2.5 higher than our "standard" predictions (cf. Fig.[3j, 
since we assume that the continuum close to Lya is unchanged 
by this departure from case B. 

A few words of caution about W(Lya) are appropriate. First, 
note that for our computations of W(Lya) we use the continuum 
flux at 1215.67 A obtained from linear interpolation of the total 
(stellar + nebular) continuum (in log) between 1 190 and 1240 A, 
chosen to avoid underlying stellar Lya absorption and other ab- 
sorption lines if present (cf. above). While W(Lya) is well de- 
fined theoretically, comparisons with observations require some 
caution, given the possible complexity of the continuous spec- 
trum shortward (due to the IGM in particular) and longward of 
Lya (due to non-monotonic shape of the nebular continuum), 
and given different choices of broadband filters (see e.g. the sim- 
ulations of Hayes & Ostlin (2006)). 

How much of the total radiative energy from the starburst is 
emitted in the Lya line? For constant star-formation the frac- 
tion of the Lya to the bolometric luminosity, L(Lya)/Lbd, is 
shown in Fig. [6] for all metallicities and IMFs. At solar metal- 
licity and for a Salpeter IMF we obtain the well-known esti- 
mate of L(Lya)/Lboi ~ 3-6 % found in the first papers pro- 
moting Lya searches at high redshift (e.g. Partridge & Peebles 
1967). The fraction of luminosity emitted in Lya increases with 
decreasing metallicity, due to the higher ionizing photon flux 
output per unit stellar mass. When case B departures at low 
metallicity are taken into account, we find that L(Lya)/Lboi can 
reach up to ~ 20-40% depending on the IMF, i.e. up to 10 
times more than expected from earlier calculations! The highest 
values are comparabl e to those from the photoionization mod- 
els of Panagia (120021) using very hot black body spectra. For 
younger populations L(Lya)/Lboi is less dependent on the IMF 
than for SFR=const shown here; values of L(Lya)/Lboi ~ 0.15- 
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0.20 (0.35-0.40) are obtained at low Z with our standard (depar- 
ture from case B) assumptions. 

3.5. He ii line emission from very metal-poor starburst 
galaxies 

Our standard predictions for He n A 1640 (using Eq. [3]l for con- 
stant star-formation and young bursts are shown in F ig. [7] They 
complement our earlier predictions in Schaerer (2003), and show 
the expected behaviour. Clearly, strong nebular He n /U640 
emission from starbursts is only expected at very low metallicity 
and for IMFs producing enough massive stars. As we will show 
below, the predicted intensity of He n /11640 (and other He + re- 
combination lines) depends, however, also on the ionization pa- 
rameter and on the ISM density to some extent. Complete pho- 
toionization models predict generally fainter He n A 1640 emis- 
sion, as discussed in Sect. [4] 
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Fig. 6. Fraction of the Lya luminosity to the total bolometric lu- 
minosity, L(Lyar)/Lboi for SFR=const as a function of metallic- 
ity and IMF. Results for different IMFs are shown using the same 
colour codes as in Fig.[T](cf. TableQ}. Thin lines show results us- 
ing our "standard" Lye predictions; thick lines the improved re- 
sults accounting for departures from case B at very low metallic- 
ity. Note the resulting strong increase of our revised L(Lya)/L\, \ 
values from solar to very low metallicity! 



3.6. Importance of the nebular continuum 

Figure |8] shows the contribution of nebular continuous emission 
to the total UV light at 1500 A for all metallicities, IMFs, and 
for the usual limiting cases of star-formation histories. Whereas 
for "normal" metallicities and IMFs the contribution is relatively 
small ($5% for SFR=const), the importance of the nebular con- 
tinuum is larger for young bursts and/ or low metallic ity, as al- 
ready stressed bv lSchaererl d2002l) and lSchaererl ((2003). 

Since the nebular continuum at 1500 A is generally dom- 
inated by the two-photon continuous emission process, depar- 



tures from Case B will lead to stronger nebular emission than 
shown here at very low metallicity (see Sect. |4). To first order, 
the total two-photon emission is then enhanced by a factor P for 
low ISM densities, increasing thus the contribution of the nebu- 
lar continuum to the total (stellar + nebular) emission. At high 
density, the two-photon emission tends to zero. However, since 
the shape of the nebular continuum depends on the detailed con- 
ditions in the nebula (density, temperature, etc.) which are not 
constant as assumed in the synthesis models, it is not possible to 
predict more accurately how departures from case B affects the 
nebular continuum, without resort to photoionization models. 

3.7. Predicted UV slope 

From our synthesis models we also measured the slope of the 
UV continuum with and without nebular emission (cf. Schaerer 
& Pello 2005). In Fig. [9] we show a condensed overview of var- 
ious /3-slopes for very young populations (ZAMS), representing 
the steepest slopes, i.e. the minimum for [}, predicted from mod- 
els. As before, the predictions are shown for all IMFs and metal- 
licities. Thick (thin) lines show ^61550 O82000X defined as the slope 
between 1300-1800 (1800-2200) A respectively The upper set 
of lines shows the j3 slopes of the total spectrum, including stel- 
lar and nebular continuum, the lower lines using the pure stellar 
spectrum. 

Clearly, the UV slope is strongly affected by nebular emis- 
sion , leading to a si gnificant flattening of the spectrum (cf. Figs. 
2 in ISchae rer 2003). While the stellar SED becomes steeper with 
decreasing metallicity, the total spectrum exhibits the opposite 
behaviour. If we assume a varying contribution of the nebular 
continuum we may obtain any intermediate value of (3 between 
the "stellar+nebular" and "pure stellar" cases. This fact and the 
depende nce of /? on the star-form ation history and age (see e.g. 
Fig. 1 in lSchaerer & Pelio1l2005l) . shows that the UV slope can- 
not be used to determine metallicity. For constant star-formation 
(not shown here), the bulk of the models show equilibrium val- 
ues of j3 ~ -2.6 to -2, quite independently of metallicity and IMF. 

Fig.|9]also shows a difference of the order of 0.2-0.3 between 
/?i50o and /?2000' Such differences may be relevant for compar- 
isons of the UV slope estimated from broad-band filters. Last, 
but not least, the precise shape of the nebular continuum is de- 
termined by the detailed nebular structure (i.e. its detailed tem- 
perature and density structure) and can be affected by departures 
from case B, as the case of the two-photon continuum discussed 
in depth below. 

4. Nebular predictions using photoionization 
models 

For low metallicity nebulae ionized by very hot stars, the conven- 
tional case B predictions for line and continuum emission are not 
good approximations to the appropriate nebular astrophysics. In 
this Section we present the predictions of detailed photoioniza- 
tion modeling for such nebulae with metallicities ranging from 
zero (primordial = Pop III) to solar, explain the origin of depar- 
tures from case B, and present a parameterisation of the results 
than can be readily employed for the interpretation of low metal- 
licity nebulae. For simplicity, we first use black body spectra as 
ionizing sources and discuss later how to compare these with 
stellar SEDs. We have examined the nebular emission as a func- 
tion of the stellar (black body) temperature T^, the hydrogen 

3 We use the standard definition F A oc ?P. 
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Fig. 5. Predicted Lye equivalent width as a function of metallicity for constant star formation (left panel) and very young (< 1-2 
Myr) bursts (right panel). Note the different vertical scales on the two plots. Thin lines show our "standard" predictions, thick lines 
the predicted W(Lya) accounting to first order for departure from case B following Eq. [8] (assuming low density, i.e. f co \\ - 2/3, 
and neglecting the increase of the two-photon continuum), leading to an increase by up to a factor ~ 1.5-2.5 at low metallicities 
(Z £ Z co u). Results for different IMFs are shown using the same colour codes as in Fig.[T](cf. Table[T]i. 
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Fig. 7. Predicted He n /il640 equivalent width as a function of metallicity for constant star formation (left panel) and very young 
(< 1-2 Myr) bursts (right panel). Note the different vertical scales on the two plots. Results for different IMFs are shown using the 
same colour codes as in Fig.Q](cf. TableQ]). Note that photoionization models predict generally fainter He n /11640 emission, hence 
lower equivalent widths, except for high ISM densities (see Sect. 14. 31 . 



number density n#, the ionization parameter (U), and the neb- 
ular metallicity (Z ne t,), using the models described above (Sect. 



Our analysis will focus mainly on the Lya line and the asso- 
ciated two-photon continuum emission, on nebular He n emis- 
sion, and on the global nebular continuum. 
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Fig. 8. Contribution of nebular continuum emission to the total emission in the restframe UV at 1500 A as derived from our 
evolutionary synthesis models. Shown are model for constant star formation (left panel) and very young (< 1-2 Myr) bursts (right 
panel). Results for different IMFs are shown using the same colour codes as in Fig.Q](cf. Table[TJ. 
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Fig. 9. Left: Predicted UV slope /61500 and for all IMFs and metallicities. The values are shown for very young (ZAMS) populations, 
which correspond to the bluest possible slopes (i.e. minimal /? values). The upper set of lines shows the UV slopes of the total 
spectrum (stellar + nebular continuum), the lower lines using the pure stellar spectrum. Results for different IMFs are shown using 
the same colour codes as in Fig.Q](cf. Table[TJ. Right: Same as the left panel, for #2000- Note the difference between fi\ 500 and #20001 
indicative of the deviation of the true SED from a simple power law. 



4.1. Lya line and two-photon continuum emission 

In Fig.[TO]we show the deviation of the Lye (black symbols) and 
2y emission (red) and other quantities as predicted from pho- 
toionization models with primordial composition and different 
densities with respect to their case B values. To calculate case B 
luminosities we adopt Eq. [2] and assume L{2y) = 0.5 c\ Q(H) 



for the luminosity in the two-photon continuum. The numerical 



factor 0.5 ~ /a^ p is appropriate for low densities. For sim- 



plicity we always compare the results to the low density limit 
case B predictions. 

At low density, the luminosity of both Lya and 2y emis- 
sion are increased by a factor of ~ 1.3 to 3 over the black body 
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Fig. 10. Lya (black), 2y (red), Ha (green), H{3 (blue), He n 1640 A (orange) luminosities over their predicted case B luminosities for 
the primordial nebula. The black solid line shows P (Eq. |6j, the dashed lines show P x / C oii//coii(log(nH = 1)) (cf. Eq. [7j. Different 
symbols represent different ionization parameters, namely open diamonds correspond to log(f/) = -1, stars to log(t/) = -2, open 
squares to log(t/) = -2, and open triangles to log(f/) = -4. It can be seen that Lya and 2y luminosities scale with P at low densities 
and Lya is further enh anced at h i gher d ensity. Case B values have been calculated using luminosity coefficients for r e =30,000 K 
and riH-10 2 cirT 3 from lSchaererl ( 2003). For simplicity we always compare the results to the low density case B limit. 



temperature range considered here. The physical reason for this 
strong departure from case B is due to collisional effects, which 
increase the population of the n — 2 level of hydrogen from 
which additional ionizations can take place, leading overall to 
an increased ionization rate in the nebula. In equilibrium, this 
implies an increased recombination rate and higher Lya- and 
2y luminosities. Collisional excitation is significant in (very) 
low metallicity nebulae with hot ionizing sources because radia- 
tive cooling is reduced, leading to higher electron temperatures. 
Photoionization from excited states (here from n - 2), in partic- 
ular, is not taken into account under case B assumptions. 



To describe quantitatively the effect of enhanced photoion- 
ization rates, including ionization from the excited n = 2 state, it 



is sufficient to consider the mean energy of the ionizing photons 
in units of the ionization potential of hydrogen, P, given by 



P = 



13.6 eV 



13.6 eV 



Fydv 



f°° 

Jl3.6 eV hv ' 



/ (13.6 eV), 



(6) 



where E is the average ionizing photon energy in the Lyman con- 
tinuum in units of eV. Indeed, as Figure [10] shows, the enhance- 
ment of the Lya- and two-photon continuum emission, which 
both measure the effective recombination rate, scales very accu- 
rately with P at low density. This scaling shows that the available 
energy of the Lyman continuum photons - which is in excess 
of the necessary minimum of 13.6 eV - is "optimally" used to 
maximise the number of photoionizations, leading to an effec- 
tive increase of the ionization rate from oc (9(H) tooefx Q(H), 
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and hence to the same increase in the recombination rates, i.e. 
also in the Lya and 2y luminosities. 

The effect of "boosted" Lya and 2y emission just discussed 
depends on the nebular metallicity. This is because more metals 
are present and more efficient cooling results in cooler gas than 
in the primordial metallicity case, decreasing the collisional ef- 
fects for H. The range over which these effects take place is dis- 
cussed below. 
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0.5 L. 



Panagia 1973 
this work 



3 4 
og(n(H) [cm"]) 



The corresponding values of /oil, obtained for our grid of mod- 
els for primordial nebulae are plotted in Figure QT] as a function 
of the input hydrogen density. Also shown is, for comparison, 
Panagia's expression (Eq.[7} assuming n p — Notice that / co n 
presented here is not derived from the Lya and 2y luminosities 
(which ranges from 2/3 to 1 as in Panagia's formula) but it is 
the fitting formula which has been calculated by comparing the 
right hand side of equation [8] with our CLOUDY results. 
For convenience, we derive the following 3-parameter fit: 



1 + a X tin 
b + c xtiH 



(9) 



with a = 1.62 X 1(T 3 , b = 1.56 and c = 1.78 x 1(T 3 . Th e coef- 
ficients obtained here differ from those of lPanagial (1 19731) due to 
the fact that we have used a suite of photoionization models cov- 
ering a wide range of and U resulting in different T e and a 
different balance between processes in the nebula. However, the 
change in / co u has little effect on the values of L(Lya) for low 
density. We use / u defined in this way for convenience since it 
allows us to reproduce the CLOUDY results. As Fig.fTTIshows. 
our fit formula allows us to describe density effects with an ac- 
curacy of +10% at high density and somewhat better at low 

The black dashed lines in Figure [10] showing P x 
/coii/^oii(log(nH = 1)), show how well our analytical predic- 
tion for Lya, including also density effects, works at primordial 
metallicity. 



Fig. 11. / co ii factor accounting for the density effects. The red 
line represents the 3-parameter fit obtained on the basis of the 
photoionization modeling (/ co ii) to be used in formula [8] to cal- 
culate L{Lya). At each density, vertical lines show the spread in 
/oil values arising from the different T/,b and U models used: 
models with higher Ty, and lower U have higher / co ii values. 
The fit has been made for primo rdial nebular abundances and 
the black dotted line is / co ii from Panagia| d!973l) . 



4.1 .1 . Density effects 

Figure [10] also shows that Lya is further enhanced at high den- 
sities at the expense of two-photon emission, the sum of the two 
luminosities being essentially constant if all other parameters are 
kept fixed. Collisional mixing of the relative populations of the 
n — 2 levels between 2 2 S and 2 2 P of hydrogen alter the rel- 
ative importance of the rates of Lya and two-photon transitions 
resulting in increased line emission relative to the 2y co ntinuum. 

To correct f or this effect u nder case B assumptions, Panagia 
(I1973L see also lStiavellil l2009) writes L{Lya) = Q(H) x hv Lya X 
/coii (cf. Eq. [JJ, where 
1 + a X n„ 



H n / H„, log(U) 



1.5 + b x n T , 



(7) 



with a — b — 1.35 x 10 4 , and where n p is the proton density. 
In this way the factor f co \\ ranges from 2/3 ~ a^ p /aB for low 
densities to /con = 1 for high densities. Under these assumptions, 
the two-photon continuum luminosity is L(2y) = (2(H) x/iVL yff x 

(Wcoll). 

To derive the corresponding numerical factor from our de- 
tailed photoionization models and to separate the density ef- 
fect from the enhancement factor P found above, we define / co n 
through: 




L(Lya) = (1 - / eS c)g(H) x hv Lya x P x /« 



(8) 



Fig. 12. The Balmer decrement as a function of T^, and n# for 
primordial metallicity case. Note that Ha is boosted for high Ty, 
due to collisional excitation of hydrogen. 



4.1.2. Applicable metallicity range 

The effect of collisional excitations, allowing for photoioniza- 
tion of H from excited states, depends on the electron tempera- 
ture of the nebula, and hence on its metallicity. To find the range 
of nebular metallicity where our analytic expressions for L(Lyo') 
holds, we have computed model grids at different metallicities. 
We can then ascertain how much our fit to the primordial case 
deviates from CLOUDY models as a function of metallicity. 

The critical metallicity where the luminosity obtained by our 
fit deviates on average (for the entire grid of models) by 10 % 
from the correct one is Z co u ~ 0.03 Z . At higher metallicities, 



12 



A. Raiter, D. Schaerer, R. Fosbury: UV properties of starburst galaxies 



cooler nebulae are produced and hence collisional effects play a 
smaller role in boosting hydrogen emission. However, the effect 
only decreases gradually for higher metallicities and depends on 
the detailed condition of the gas. For instance, in our simulations 
for Z ne b = 0.05 Z , the average deviation of the results obtained 
with our formula (primordial case) with respect to the CLOUDY 
results is 13%, while for Z ne b = 0.1 Z Q it is 20%. The value of 
the "transition" metallicity Z co n ~ 0.03 Z below which Lya 
(and the 2y continuum) are significantly boosted (modified) by 
departures from case B, should therefore only be taken as an 
indicative value. Tailored photoionization models are necessary 
for more accurate predictions. 



4.2. Other H lines 

The luminosities of other hydrogen lines are close to their case B 
predictions as is shown in Figure[10]for Ha (green symbols) and 
H/3 (blue symbols). They are basically independent of U. 

These lines are used to calculate the Balmer decrement 
which is conventionally used to measure the extinction. The 
Balmer decrement (Ha/H/?) obtained from our grid of models 
for log(f/) = -1 is shown in Figure [12] It increases towards 
higher 7bb and exceeds the values commonly used, for ex- 
ample 2.86 for n e = 10 2 cm 3 and T e = 10 kK (Dopita & 
Sutherland 1996). The reason for the enhanced HafH/3 ratio is 
the collisional excitation of hydrogen, one of the contributors 
to enhanced Lya and 2y emission. This effect (see Osterbrock 
2006) occurs when photoelectrons carrying enough energy when 
colliding with H atoms, excite higher levels followed by radia- 
tive cascade. Enhanced Ha emission increases the Balmer decre- 
ment which must be accounted for when calculating the extinc- 
tion. An enhanced Balmer decrement was already found, e.g. in 
tailor ed photoionization models for the metal-poor galaxy I Zw 
18 dDavidson & Kinmarifl^lStasinska & Schaerer|[T999l) and 
other giant H n re gions (cf. iLuridiana et al.ll2003h . and has also 
been addressed bv lLuridianal d2009» ! 



4.3. He ii lines 

The He n 1640 A line is of a particular interest since it has been 
considered as one of the signatures of Pop III or very metal- 
poor stars, which are expected to have a very hard ionizing spec- 
trum emitting copious He + ionizing photons. Figure [TUl (orange 
symbols) shows He n 1640 A emission line luminosities rela- 
tive to their case B values, calculated for the same T e and n e 
as the hydrogen lines. For lower ionization parameters, the line 
becomes weaker. Th i s is du e to an effect already discussed by 
Stasihska & Tvlendal (fl9 86) in the context of planetary nebulae. 
Photons with sufficient energy to ionize He + are also able to ion- 
ize hydrogen. When calculating the luminosity of He n in synthe- 
sis models, it is assumed that every photon with an energy >54 
eV ionizes one He + ion. The absorption of some high energy 
photons by H atoms results in a decrease in the number of pho- 
tons available for producing He n emission. The effect becomes 
significant at low ionization parameters since then the He + /H° 
fraction decreases and the probability of the high energy photon 
being absorbed by a H atom rather than by He + is higher. This 
means that the synthesis models give an upper limit for the He n 
luminosity. We have not found a simple analytical prescription 
to account for the effect of this process on the He n luminosity. 

The behaviour of the equivalent width of He n is discussed 
in Sect. [5] 
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Fig. 13. Continuum of T^= 100,000 K models as a function 
of the nebular metallicity for a constant density and ionization 
parameter. The red solid line represents the transmitted stellar 
(black body) continuum. The thin lines are nebular only contri- 
butions to the continuum. The thick lines are the total contin- 
uum emission (stellar + nebular). The first model is normalized 
at 1240 A and the rest are rescaled accordingly. 
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Fig. 14. Continuum of T^= 100,000 K models as a function of 
the hydrogen number density at the primordial metallicity and 
constant ionization parameter. The red solid line represents the 
transmitted stellar (black body) continuum. The thin lines are 
only nebular contributions to the continuum. The thick lines are 
the total continuum emission (stellar + nebular). At higher den- 
sities 2y continuum gets destroyed due to collisions. First model 
normalized at 1240 A and the rest rescaled accordingly. 



4.4. Nebular continuum emission 

Nebular emission arising in H n regions can contribute signif- 
icantly to the measured spectrum. The total nebular continuum 
is the sum of the free-free, free-bound and 2y continua of hy- 
drogen and helium. Several nebular parameters affect its shape. 
It depends on the electron temperature, ionization parameter, 
nebular metallicity and particle density in the gas, as e.g. illus- 
trated by Bottorff et al. (2006). In Figure [14] we show the den- 
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sity dependence of the nebular continuum, which mainly affects 
the 2y emission and hence also the shape of the continuum be- 
tween Lya and the Balmer jump. Figure [13] shows the depen- 
dence on the nebular metallicity. For higher metallicities, the in- 
creased cooling lowers the total emission in the nebular contin- 
uum, and its shape is altered as expected from atomic physics 
(see Bottorff et al. 2006) due to the decrease of the average 
electron temperature at higher metallicity. Evolutionary synthe- 
sis models such as ours, cannot properly describe the variety of 
shapes and strengths of the nebular continuum shown here, since 
they rely on simple assumptions such as constant nebular den- 
sity and temperature. Again, tailored photoionization models are 
necessary for more accurate predictions. 

5. Photoionization models for realistic SEDs 

In this Section we show how to relate realistic SEDs obtained 
e.g. from evolutionary synthesis models with the results from 
the photoionization models discussed above, which used black 
body spectra. Finally, we show the updated predictions for the 
equivalent widths of Lya and He n /11640 lines obtained from 
our CLOUDY models. 

5. 1. How to connect realistic SEDs with black body 
calculations 

Since the number of H ionizing photons and their mean energy 
are the main parameters determining the luminosity of the H re- 
combination lines and of the 2y continuum, Q(H) and P suffice 
to compute these quantities from arbitrary SEDs. To account also 
for density effects, Eqs.[8]and|9]can be used to yield the correct 
Lya luminosity. Numerical tests using SEDs described in Sect. [3] 
show the validity of this result, typically within 2-3 %. For black 
bodies the relation between P and r bb is given by 

r bb [kK] = -53.42 + 69.85 P (10) 

to a good approximation. The black body temperature r bb (P) 
corresponding to the more realistic SED is thus easily deter- 
mined. 

For other quantities, such as the He n line luminosity, the 
correspondence between black bodies and other SEDs is differ- 
ent. For example, since the relative line ratios of He n/H lines 
depend to first order on the relative number of He + /H ionizing 
photons, the black body with the same hardness Q(He + )/Q(H) is 
the most appropriate. For black body spectra one has, to a good 
approximation: 

r bb [kK] = 314.6 + 382.4x + 268.7 x 2 + 

+ 103.3 Jt 3 + 20.12 x 4 + 1.545 x 5 . (11) 

where x - log(£)(He + )/(2(H)). Again, we have tested a number 
of models and confirmed that one can get the same luminos- 
ity (within a few %) using SEDs and corresponding black body 
models, for otherwise identical nebular parameters. Differences 
of ~ 10% can appear in case of significantly (several eV) dif- 
ferent values of P. In principle not only the stellar Q(He + )/Q(H) 
ratio determines the He n emission, but also the conditions in 
the gas. Since P establishes the electron temperature, significant 
differences in P result in different luminosity coefficients (re- 
combination rates and emissivities). 

In Fig. [15] we show such a comparison between a realistic 
SED from our synthesis models (Pop III, Salpeter IMF 1-100 
M at zero age, instantaneous burst) and the appropriate black 
body, chosen such as to reproduce the correct H line luminosities 



and nebular continuum. Here the average photon energy in the 
Lyman continuum is E — 30.033 eV, hence T bb = 102 kK. The 
input spectra, shown on the left panel, are scaled to the same 
total Lyman continuum flux Q(H). The nebular parameters are 
log(«n) = 1 crrT 3 , log({/) = -1, and primordial composition. 

Several interesting points are illustrated with this figure. 
First, as the left panel shows, it must be remembered that realis- 
tic SEDs show a Lyman break due to the hydrogen opacity in the 
atmosphere of the hot stars responsible for the flux. For young 
ages this break corresponds typically to a flux increase by ~ 0.2- 
0.3 dex at 912 A, as discussed e.g. by Schaerer (20030 In con- 
sequence, calculations relying on black body spectra underesti- 
mate the observable UV continuum, as shown in the right panel 
of Fig. [15] leading to non-negligible differences e.g. for equiva- 
lent widths predictions of emission lines (cf. below). In Fig. Q3] 
(right panel) we also show the stellar + nebular SED predicted 
by our synthesis models assuming constant nebular density and 
temperature and case B (green line). The differences between 
this and the full CLOUDY model are mainly due to the boosting 
of the 2y continuum discussed earlier, and to varying n e and T e . 

5.2. Lya equivalent width predictions 

The measurement of the equivalent width of Lya is both obser- 
vationally and interpretationally difficult, since the continuum 
around it is affected on one side at high redshifts by Lya for- 
est absorption and on the other by an unknown combination of 
2y continuum and starlight. Furthermore, when not measured 
from spectroscopy, it is common practice to measure W(Lya) 
via a line to continuum ratio with the continuum estimated at 
another (usually longer) wavelength Acorn e -g- Wa = i-iine/^w- 
Different methods to m easure W(Lya) have been reviewed by 
lHaves & Ostlinl d2006l) . To illustrate the possible impact or un- 
certainty related to the way the continuum is estimated we subse- 
quently plot Lya equivalent widths predicted from the CLOUDY 
models using the continuum at both 1240 and 1500 A as a refer- 
ence. 

Figure [16] shows the predictions for W(Lya) from all 
CLOUDY models for primordial composition, as a function of 
r bb , ionization parameter and hydrogen density. Black sym- 
bols show the predictions using black body spectra and account- 
ing for the nebular continuum; red symbols the same but ne- 
glecting the nebular continuum. Blue symbols show the results 
from CLOUDY models using SEDs from our Pop III synthesis 
models for zero age populations, plotted at the corresponding 
r bb (P) value (see above). Pink symbols show W(Lya) predicted 
from our standard synthesis models neglecting the boost of Lya 
(case B departure). 

Overall W{Lya) increases with the r bb or equivalent since 
the Lya luminosity (oc Q(H)xP) increases more rapidly than the 
continuum flux close to Lya. For very hot models and if 1500 
A is taken as a reference for the continuum, W(Lya) tends to 
a maximum value (here ~ 2000 A, as shown in col. 2) since 
line emission and the dominating nebular continuum scale in the 
same manner. For the reasons discussed above (absence of the 
Lyman break), predictions based on black body spectra (black 
and red symbols) overestimate W(Lya) compared to models in- 

4 Th e amplitudes of the Lyman break, 912 + /912~, given in Schaerer] 
(2003) for very low metallicities are correct, but misleading. Indeed, 
with the definition adopted there, the 912 + flux average (over 1080- 
1200 A) includes strong absorption from the Hen /11084 line. This 
leads to values of 912 + /912~ < 1 despite the fact that all models show a 
Lyman break in absorption. 
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Fig. 15. Left: Comparison of input stellar spectra for photoionization models showing a realistic SED (blue line; Pop III, Salpeter 
IMF 1-100 M at zero age) and the corresponding black body spectrum (red line) with the same average photon energy {E = 30.033 
eV) in the Lyman continuum. The SED model is normalized at 1240 A and the rest rescaled to match the same ionizing flux Q(H). 
Right: Input stellar spectra (dash-dotted, same colours as left panel) and predicted nebular and total continua from our CLOUDY 
models for log(n(H)) = 1 cm~ 3 and log(t/) = -1 (solid lines, using the same colours as in left panel; cf. also inset for symbols). The 
green solid line shows the total (stellar + nebular) SED from our evolutionary synthesis model assuming constant nebular density 
and temperature and case B. Notice the difference between case B and the CLOUDY result due to enhanced 2y emission. 



eluding more realistic stellar SEDs (in blue). Differences be- 
tween black body models including or neglecting the nebular 
continuum (black and red) decrease with increasing ISM den- 
sity, because of the decrease of the 2y continuum. This also 
explains the decreasing difference between W(Lyo') using 1240 
and 1500 A as a reference, when increases. Last, but not 
least, W(Lya) predicted by our standard synthesis models for the 
ZAMS (pink symbols) fall below the more accurate CLOUDY 
predictions (blue), since the latter account for case B departures. 
Once applying the simple correction given by Eq.|4]one obtains 
the result which takes into account the increase in Lye line flux. 
However, that is only the first order correction which does not in- 
clude the difference coming from the continuum, mainly a sim- 
ilar increase in the 2y continuum. In cases when the continuum 
is measured around Lya (1240 A) stellar flux is dominant so the 
correction described above gives a good estimate of the correct 
W(Lya-) value (see Fig.fTTt. 

5.3. He ii A1 640 equivalent width 

The right column in Fig. [16] shows the predicted He n /U640 
equivalent widths for the same CLOUDY models except that 
Tbb(^) from Eq. QT|is used here . The nebular continuum pre- 



dicted at 1640 A is included in all calculations of W(1640). 
In contrast to Lya, the He n equivalent width depends on 
the ionization parameter, since the line luminosity changes 
due to the "Stasinska-Tylenda effect". For high ISM densities 
W(1640) increases due to the decrease of the nebular contin- 
uum. Differences between models assuming black body spectra 
(black symbols) and realistic SEDs (blue) are minor for W(1640) 
as long the nebular continuum dominates at these wavelengths. 
This is true for sufficiently hot spectra (> 80 kK). Compared 
to predictions from our standard evolutionary synthesis mod- 
els (pink symbols) the He n Al 640 equivalent widths predicted 
by CLOUDY are significantly lower for T hb > 100 kK. This is 
due to two effects, the Stasinska-Tylenda effect reducing the line 
flux, and the increase of the 2y nebular continuum. At lower 7bb 
, case B predictions give an upper limit (corresponding to our 
highest ionization parameter case) for He n /U640. Only for a 
high ionization parameter and/or high density do photoioniza- 
tion models pre dict W(1640) values as high as those given by 
ISchaereri(l200l . 
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Fig. 16. Equivalent widths (W,\) of Lya at 1240 (column 1) and 1500 A (column 2) and He n /11640 (column 3) for a primordial 
metallicity nebula at different densities and ionization parameters (marked with different symbols as defined in the inset). The red 
symbols show an upper limit for W^Lycr) since they are calculated using only stellar (black body) continuum. Blue symbols - the 
most realistic predictions - show the results obtained for the corresponding stellar population SEDs (models with the same P in the 
case of Lye and with the same log((2(He + )/2(H)) for the He n analysis, for simplicity only models with log(t/) = -1 are shown) 
and the pink ones their case B predictions. Notice the dependence of the He n /U640 equivalent width on the ionization parameter 
and the discrepancy between its case B and photoionization model predictions for real SEDs. 



5.4. Summary: correction of evolutionary synthesis model 
results for departures from case B 

In this section we give a recommended step-by-step procedure 
which allows the correction of the synthesis models for the de- 
partures from case B in ionization bounded regions (i.e. with no 
leakage of Lyman continuum photons), relevant especially for 
the case of very hot stars and very low metallicities. 
For L(Lya): 

- Take the SED (function of age, IMF, ...). 

- Get (or compute) the mean energy of Lyman continuum pho- 
tons, P, in units of 13.6 eV. 

- Compute the "collisional" factor ^ n from our fit formula 
(Eq.lU. 



— Multiply the synthesis model result L(Lyff) by P */ co ii/(2/3) 
(Eq.gll. 

For W(Lya-): divide the equivalent width obtained by the syn- 
thesis model of interest by 2/3 and multiply it by P * / co u. This 
provides the first order correction for W(Lya), neglecting the in- 
crease of the nearby nebular continuum. An example following 
this procedure is shown by the green symbols in Fig.fTTl 

For He n recombination lines: no simple correction is possi- 
ble. Synthesis models give an upper limit for the luminosity, an 

upperlimit for equivalent width up to log(Q(He + )/2(H)) 1.2 

for the primordial metallicity case and overestimate it for the 
highest Q(He + )/Q(H) cases (low density) due to the lack of ac- 
counting for enhanced 2y emission. 
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Fig. 17. Equivalent widths (W A ) of Lya at 1240 A (first col- 
umn, second row of Fig. [ToT l. All the symbols as in Fig. [16] 
Additionally, the green symbols show the case B SED results 
(pink symbols) corrected using our formula (first order correc- 
tion; only the enhancement in the line flux is taken into account, 
there is no correction for the continuum). The difference coming 
from not correcting for the nebular continuum enhancement is 
visible for the hottest model. 



6. Discussion and implications 

6.1. Dependence on model assumptions 

All our CLOUDY models have been computed for ionization 
bounded nebulae corresponding to no escape of Lyman con- 
tinuum photons, i.e. f esc = 0. For density (or matter) bounded 
cases with / esc > 0, quantities such as hydrogen recombina- 
tion line luminosities and the luminosity of the nebular contin- 
uum longward of Lya- can be rescaled to first order by scaling 
with (1 - /esc), as already introduced in Eq. [8] Such a scal- 
ing is appropriate since the corresponding emissivities are ap- 
proximately constant across the nebula. However, since He + 
recombination lines (if present) are emitted only in the inner- 
most parts of the nebula, He n line luminosities will not be re- 
duced in density bounded regions, as long the radius of the He + 
sphere remains smaller than that of density bounded H n region. 
These behaviours ar e also expected in no n-spherical regions, as 
demonstrated e.g. bv lJohnson et al.ld2009l) . In consequence, both 
L(Lya) and W(Lya) will decrease with increasing Lyman con- 
tinuum escape (although the latter not proportionally with / e scX 
whereas the opposite will be true for W(He n A 1640) and for 
other He n lines, where the underlying continuum is mostly of 
nebular origin. In other words, the decrease of W(He n A 1640) 
found in this paper due to competition with H ionizing photons, 
may be mitigated in objects with significant leakage of Lyman 
continuum radiation, lead ing ag ain to higher He n equivalent 
widths (cf . ISchaererll2002l |2003b . For other effects of "leaking", 
i.e. star-forming galaxies with density boun ded H n regions, see 
the photoionization models of llnoud(l2010l) . 

We shall now discuss some implications of our modeling for 
the interpretation of the hydrogen and helium emission lines and 
the nebular continuum observed in high redshift emission line 
galaxies. 



6.2. Lyman continuum output 

As metallicity decreases, the Lyman continuum output of the 
stellar population, both relative to the observed stellar UV con- 
tinuum longward of Lya and per unit stellar mass, increases. The 
amount of this increase depends obviously also on the detailed 
shape of the IMF. For example, for a Salpeter IMF with a uni- 
form upper mass cut-off of 100 M Q , the ionizing photon rate per 
unit UV luminosity, Qh/£i5oo> increases by a factor 3 approxi- 
mately between solar an zero metallicity. For the most extreme 
IMF considered here this increase is up to a factor 10 (see Fig. 
[TJ. The effect of this is to increase the line and nebular con- 
tinuum emission relativ e to any ob servable stellar continuum. 
In extreme cases (e.g., iRaiter et all 1201 Oh the nebular emission 
may completely dominate the UV/visible/NIR spectrum. Also, 
if low metallicity objects and/or different IMFs are relevant for 
the sources of cosmic reionization, their intrinsic Lyman contin- 
uum flux, generally estimated from UV restframe observations, 
may need to be revised accordingly. 



6.3. SFR(UV) 

The estimation of SFR from the level of the UV continuum (cf. 
lKennicutlH l998) is based on evolutionary synthesis models as- 
suming "standard" stellar populations and IMFs. At metallici- 
ties below solar, the output of observable UV continuum radi- 
ation (i.e. typically at 1500 A restframe) per unit stellar mass 
increases somewhat, due to the increasing average temperature 
of the stars (cf. Fig. |4j). A maximum in L1500 is typically reached 
at Z/zQ ~ 1/100. Below this, the average spectrum shifts more 
strongly into the Lyman continuum, leading to a decrease of stel- 
lar radiation in the observable UV, which is, however, compen- 
sated by increasing nebular emission. The net result, illustrated 
in Fig. [4] for constant SFR, is that star-formation rates derived 
from the UV may need to be revised downward at l ow metal- 
licities with respect to the "standard" calibration of Ken nicuttl 
( 1998), but typically by less than a factor 2. 

6.4. Lya 

One of the main findings of our study is the enhancement of 
Lya emission at low metallicities with respect to the commonly 
used case B value. This departure from case B occurs even at 
moderaly low metallicity, say Z/Z 5= 1/10, and becomes more 
important for lower Z. In consequence, higher intrinsic Lya line 
luminosities relative to the UV continuum (to Lboi, SFR, or to 
another measure of the rate of massive star-formation) are ex- 
pected and higher Lya- equivalent widths. Even higher values 
of W(Lya) can be obtained in regions with a high ISM density, 
where Lya emission can further be "boosted" at the expense of 
the 2-photon nebular continuum (see Fig. [TOt . 

With respect to our earlier predictions (cf. Sch aererll2003l) . 
W(Lya) is increased by a factor 2-2.5 for Pop III objects, and 
by more than 70% at moderately low metallicity (Z/Z Q ^ 1 /20, 
cf. Fig. [5}, for a given IMF. This increase of Lya may help to 
understand objects with large Lya equivalent widths found in 
some surveys, without the need for recourse to unusual IMFs or 
to a clumpy ISM (cf. iMalhotra & Rhoadsll2002T; iDawson et al.l 
l2004tlFinkelstein et al.ll2008l) . 

Finally, the enhanced Lya emission found here should repre- 
sent good news for searches for very distant/early galaxies, since 
the intrinsic Lya emission of metal-poor objects is shown to be 
considerably higher than previously thought, reaching up to 20- 
40% of Lboi in Pop III dominated objects. Of course, the intrinsic 
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Lya emission (and hence also W(Lya)) will be lower in objects, 
having a non-zero escape of flux from the Lyman continuum. 
Furthermore, several processes exist (e.g. dust absorption, scat- 
tering out of the line-of-sight in the ISM and/or IGM), which 
will reduce the Lya emission on the way from the source to the 
observer. 



6.5. He ii emission 

The reduced He n emission found in this paper, compared to 
our earlier predictions assuming simple case B recombination, 
may indicate that ongoing searches for signatures of Po p III stars 
using this feature (see e.g. review bv lSchaererl l2008) could be 
more difficult than anticipated so far. The corresponding upper 
limits of P op III star-formatio n rate density derived e.g. from the 
survey of Nagao et al. (2008) need then to be revised upwards. 

The detection of the He n A 1640 line indicates a significant 
number of high energy photons in the ionizing spectrum. This 
line has therefore been sought, but not detected so far at high-z 
(Nagao 2008). Prescott et al. (2009) have detected this line from 
a z=1.67 spatially extended nebula (or Lya blob) but have not 
been able to determine the nature of the ionizing source. We find 
that the line (and therefore the equivalent width) could be fainter 
than expected previously, due to its dependence on the ionization 
parameter. This suggests that deeper observations will be needed 
to detect it. Additionally, another difficulty comes from the fact 
that for the hottest (the highest g(He + )/g(H)) ratio) stars the 
strongest enhancement in 2y continuum is predicted which fur- 
ther decreases the measured W(1640). However, if Pop III stars 
are present in objects with significant escape fractions in the 
Lyman continuum, the He n A 1640 equivalent width may again 
be stronger, due to the reduced nebular continuum (cf. Sect. 16. II ). 



6.6. Nebular continuous emission 

The overall nebular continuum emission is dependent on several 
nebular parameters and can result in a variety of spectral shapes 
(Fig. [13] and 11411 . Note that our simulations have been carried 
out in the absence of dust which would add another parameter 
affecting the shape of the continuum. 

Proper photoionization modeling shows the importance of 
the 2y continuum produced in a nebula ionized by very hot 
stars. At low ISM density and primordial/low metallicity, it is 
boosted in the same way as Lya and can completely dominate 
the nebular emission at 1216-1600 A. This can affect the broad- 
band flux measurements and reduce the equivalent widths of UV 
emission lines, particularly He n. The strength of the nebular 
continuum and the expected variations of its spectral shape also 
indicate that measurements of the Lya equivalent width from 
photometry may be m ore uncertain than naively expected (cf. 
lHaves & Ostlinll2006l) . 

As can be seen from Fig. [15] the slope of the UV contin- 
uum between ~ 1300 - 2000 A obtained from our CLOUDY 
model is very similar to that predicted by the evolutionary syn- 
thesis models assuming the same case B assumption, constant 
electron density and temperature. When the 2y continuum is in- 
cluded, one cannot obtain significantly steeper slopes than found 
for this Pop III simulation, since the el ectron temperature w hich 
affects the shape of 2y emissivity (cf. iBottorff et al.ll2006h be- 
comes shallower as metallicity increases. 



7. Conclusions 

Building on the earlier calculations o flSchaererl([2002ll200l . we 
have computed evolutionary synthesis models for a wide range 
of metallicities from Pop III (zero metallicity) to solar metallic- 
ity covering a wider range of IMFs - including power laws and 
log-normal IMFs with different characteristic masses - than be- 
fore. Using these synthesis models, we present the expected UV 
properties of star-forming galaxies, including their Lyman con- 
tinuum fluxes, UV luminosity, and properties of the continuum 
(e.g. /J-slopes), as well as properties of important emission lines 
such as Lya and He n /11640 (see Sect. [3j. 

To investigate possible departures from the simple case B 
recombination theory assumed in many synthesis models, in- 
cluding ours, we have computed large grids of CLOUDY pho- 
toionization models for zero and low metallicity nebulae, using 
both black-body spectra and the SEDs predicted by our synthesis 
models (Sects. [4]and[5]l. 

Our main conclusions from the photoionization models are 
the following: 

- Due to departures from case B (collisional excitation and 
ionization from excited levels), both Lya and 2y continuum 
emission can be significantly enhanced at low nebular densi- 
ties. Their strengths are found to scale with the mean photon 
energy of the ionization source in the Lyman continuum. 

- The equivalent width of Lya can be larger than expected 
from case B calculations due to the line flux enhancement. 
The measurement can also be affected by 2y emission if the 
continuum is measured at longer wavelength. 

- Further enhancement of Lya at the expense of 2y emission 
can occur due to collisional mixing between the hydrogen 
2 2 S-2 2 P levels at higher densities. 

- He n emission line fluxes (and consequently their equivalent 
widths) can be significantly decreased due to their depen- 
dence on the ionization parameter. This could make searches 
for the He n A 1640 line at high-z more difficult. 

- The enhancement of the 2y continuum and its dependence on 
nebular parameters can result in reduced equivalent widths 
of the UV emission lines and also change the UV restframe 
colours of high-z galaxies. 

Our results are of relevance to searches for, and interpreta- 
tion of observations of, metal-poor and/or high-z galaxies which 
may host the first stars that appeared in the Universe. 
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